Nonlinear microrheology of dense colloidal suspensions: a mode-coupling theory 
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A mode-coupling theory for the motion of a strongly forced probe particle in a dense colloidal 
suspension is presented. Starting point is the Smoluchowski equation for N bath and a single probe 
particle. The probe performs Brownian motion under the influence of a strong constant and uniform 
external force F cx . It is immersed in a dense homogeneous bath of (different) particles also perform- 
ing Brownian motion. Fluid and glass states are considered; solvent flow effects are neglected. Based 
on a formally exact generalized Green-Kubo relation, mode coupling approximations are performed 
and an integration through transients approach applied. A first-principles theory for the nonlinear 
velocity-force relations of the probe particle in a dense fluid and for the (de-) localized probe in 
a glass is obtained. It extends the mode coupling theory of the glass transition to strongly forced 
tracer motion and describes active microrheology experiments. A force threshold is identified which 
needs to be overcome to pull the probe particle free in a glass. For the model of hard sphere parti- 
cles, the microscopic equations for the threshold force and the probability density of the localized 
probe are solved numerically. Neglecting the spatial structure of the theory, a schematic model is 
derived which contains two types of bifurcation, the glass transition and the force-induced der- 
ealization, and which allows for analytical and numerical solutions. We discuss its phase diagram, 
forcing effects on the time-dependent correlation functions, and the friction increment. The model 
was successfully applied to simulations and experiments on colloidal hard sphere systems [I. Gazuz 
et. al., Phys. Rev. Lett. 102, 248302 (2009)], while we provide detailed information on its derivation 
and general properties. 



PACS numbers: 

I. INTRODUCTION 

Tracking the diffusive motion of a colloidal tracer par- 
ticle embedded in a soft material can be used to test its 
linear rheological [T] properties. Corresponding experi- 
mental techniques were developed during the last years 
[2H5] and are now widely known as microrheology [6] [7] . 
These techniques utilize the fluctuation-dissipation the- 
orem |8), which connects the linear response of an ob- 
servable to external fields with the corresponding time- 
dependent equilibrium correlation function. 

In order to observe nonlinear phenomena in a mi- 
crorheological experiment, the tracer has to be actively 
pulled by means of an external force. Such experiments 
were performed using magnetic forces [9] [10] as well as 
optical tweezers [ITHI5] . 

The nonlinear microrheology of dense polydisperse sus- 
pensions, which are able to undergo a glass transition 
appears to be a quite interesting tool to study the lo- 
cal properties of glassy liquids and to further reveal the 
microscopic mechanisms of vitrification [9]. In Ref. [14], 
the main results of a theory of active microrheology in 
dense colloidal suspensions close to the glass transition 
were announced and compared to simulations and exper- 
iments. The present paper is intended to give a detailed 
presentation of the theory. 

We adopt and develop the ideas of the mode-coupling 
theory (MCT) [TS] HB] approach as recently general- 
ized to macrorheology [T7lll9| . In these works, a simple 
shear flow was imposed to the suspension. The theory 
succeeded in explaining such typically nonlinear effects 



as shear-thinning and yielding in a unified microscopic 
framework. 

In the case of microrheology, the external force field 
is applied to a single colloidal tracer particle. Our the- 
ory predicts the effects analogous to the yielding and the 
shear-thinning effects in macrorheology, even though the 
technical details turn out to be quite different, as com- 
pared to macrorheology. These effects include the de- 
pinning of the tracer from the surrounding glassy ma- 
trix (so that the cage of neighbours trapping the par- 
ticle gets broken) if the applied force exceeds a certain 
threshold value. Another effect is the decay of the tracer 
friction coefficient with increasing external force. The 
force-dependent probability density of a localized probe 
can be studied in detail within the present approach, as 
was done recently for low forces [20]. We extend these 
studies to higher forces. We argue that it is possible 
to simplify the microscopic approach by neglecting the 
spatial structure. This yields schematic models which 
can readily be solved numerically and analytically. The 
nonlinear probe velocity-force relations of the most sim- 
ple schematic model, which we derive in detail here and 
where only fluctuations in the direction of the force are 
included, were compared to experiments and simulations 
in [14J . Including fluctuations perpendicular to the forc- 
ing directions, the schematic model was extended in |21) . 
and discussed in detail in [23] The latter model also could 
be extended to predict force-induced diffusion parallel 
and perpendicular to the external force [53], based on 
the microscopic memory kernels which we derive here. 

This paper is organized as follows. Sec. [TT] intro- 
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FIG. 1: The colloidal "tracer" particle is pulled through the 
suspension of "bath" particles by means of the external force 
F em . The tracer radius is a s , and a bath particle has radius 
a. 



duces the formalism of our theory. The generalized 
Green-Kubo relation is derived, valid for the nonlinear 
response. We derive the MCT equations for the time- 
dependent tracer density mode correlators, which appear 
in the mode-coupling approximations for the response 
functions. In Sec. MB we study the long time limit 



of the tracer density correlator as well as its Fourier 
back transform, which can be related to the shape of 
the "cage" of nearest neighbours surrounding the tracer. 



Sec. IV is devoted to the simplified schematic models for 
the tracer density correlator. We discuss the correlators 
of the schematic models, the corresponding scaling laws 
and the resulting friction coefficient in detail. 



II. THEORY 

A. Basic microscopic equations 

The Smoluchowski equation will provide the basis for 
all the considerations in this article: 



(i) 



where fi is the Smoluchowski operator. Eq. ([fj describes 
the time evolution of the (N + l)-particlc configuration 
space probability density ^(ri, . . . , rjv, r a , t) on a coarse- 
grained time scale, i.e. it is assumed that the velocity 
fluctuations relax much faster than the configurations. 
The particles are colloids performing Brownian motion 
with diffusion coefficients Di in a Newtonian solvent. The 
particle diffusion coefficients obey the the Stokes- Einstein 
relation 



A 



67t?7 a, ' 



(2) 



where rj is the solvent viscosity and the radius of par- 
ticle i. The colloids are allowed to interact by means of 
the potential forces Pj = — diV{r\, . . . , rjy) (di denotes 
the partial derivative <9/<9i"j), whereas the hydrodynamic 
interactions will be neglected. 

We consider a single, distinguished particle (the 
"tracer" ) with position r s and the diffusion coefficient D s 



surrounded by N identical particles (the "bath"), which 
have the diffusion coefficient Dq. The tracer is pulled by 
means of the external force F ex (see Fig. [TJ through the 
suspension. The full Smoluchowski operator 

n = n + An (3) 

consists of the unperturbed part 

n = D 9 i -(d i --^-F i )+D s d a -(d s --^—F 3 ) 



i=l N 



and the perturbation due to F c 



Afl = - 



F 8 



(4) 



(5) 



F ex will be assumed to be constant in space and time. 
From now on, we set k B T = I in the Smoluchowski op- 
erator to simplify the notation. For comparisons with 
experiments or simulations, the factor k B T will be rein- 
troduced. 

In the following, equilibrium-weighted averages 
J dT^ eq . . . will appear, which will be denoted by (...), 
where 



1 



-V({r t },r s ) 



(6) 



is the equilibrium distribution of the unperturbed system 
with the statistical sum Z. We also introduce the usual 
equilibrium- weighted scalar product, which is defined as 

(A\B) = J aT^ eq A*{T)B(T) (7) 

for two configuration-space observables A and B. 



B. Nonlinear response to the external force on the 
tracer 

Let us consider the following situation. For times 
t < 0, the system is equilibrated and there are no ex- 
ternal fields. At t = 0, the external force on the tracer 
is switched on, driving the system out of equilibrium. 
Instead of assuming that the perturbation is small, like 
it is done in the linear response theory, we consider the 
general case of arbitrarily large forces. 

With the initial condition 



*(t = 0) = * eg , 



(8) 



the solution of the Smoluchowski equation can be written 
down as 



Using the operator identity 



e nt = I + I dt'e ut n 



(9) 



(10) 
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and noting that il^ eq — A£Vi> eq , we get 



#(i) = ^ eq + [ dt'e nt 'ACl^ e 
Jo 



(11) 



The mean value of an observable A(r\, . . . , r/v) at time 
t is given by: 



(A(t))= / dr*(r s t)A(r), 



(12) 



where T is a phase space point and the integration goes 
over the entire phase space. Using eq. (Ill, we obtain 



(A(t)) = (A) - D s J dT A(T) J dt'e nt '(F ex -d s )* eq . 

(13) 

We note that d s ^ eq — F s ^ eg , introduce 

fit = Cll + AO f , (14) 

the adjoint of Q with respect to the unweighted scalar 
product [21], with 

N 

fij = £ £>o (di + Pi) • di + D s (d s + F s ) • d s (15) 



and Af2^ = D s F ex ■ d s , and finally arrive at 

(A) (t) = (A) - D s F ex ■ f di! ( F s e ot l ' A ] 
Jo 



(16) 



Eq. ( 16 1 represents the generalized nonlinear Green- 



Kubo relation for the response of an observable A to the 
perturbation by the external force on the tracer. In con- 
trast to the well-known linear response expression, the 
full Smoluchowski operator f2 containing the external 
force, instead of just the unperturbed one enters eq. ( 16 1. 



The presence of external force renders the opera- 
tor ( |14[ ) nonhermitian with respect to the equilibrium- 
weighted scalar product Q. Its adjoint is given not by 
SV but by 



n°* = nl-D s F ex ■ (f s + a.) 

(the calculation is presented in the Appendix [A| 

C. Tracer mobility 

The (long-time) tracer mobility is defined as 

(v.)(t) 



(17) 



/! = lim 



t-J-oo F P 



(18) 



where v s is the tracer velocity. In the framework of 
the Smoluchowski dynamics, where the particle motion 
is overdamped, the tracer velocity v s = d t r s — &r s is a 
function on the configurational space and is given by 



v s = n s {F ex +P a ), 



(19) 



where fj, s = D s /(fcsT) is the single-particle tracer mobil- 
ity. Since F ex is given externally and has no dependence 
on the phase space of the system, the problem reduces 
to calculating the average of the force F s from the bath 
particles on the tracer. To determine (F"), the cc-th com- 
ponent of the vector (F s ), we use eq. (16), note that the 
equilibrium average (F s ) vanishes and obtain 



dt' (F. 



e nU 'F?) 



(20) 



Let us introduce the coordinate system such that F ox 
points in the positive z-direction. Then we have F GX • 
F s = F ex Fl in eq. (20). Expression (20) can be simpli- 



fied further if we employ the rotation symmetry around 
the z-axis, which our system obviously exhibits. After 
such a rotation, the phase space integral in |20] ) should 
remain the same. On the other side, rotations by angle 
7r change the sign of both the x- and the y-component 
of F s . This means that the correlators (F s z e n+ F£ ) and 
(F*e nU 'Fy) vanish and wc arc left with 



-D s F ex f dt 1 (F s z e nU ' F; 
Jo 



(21) 



As was anticipated in eq. (18), the mean force exerted 



from the bath on the tracer is parallel to the external 
force F ex . 



Expression (21) includes the force- force correlator 
C(t) = (F z e n t Ff) and leads for the tracer mobility 
to the result 



£« = Ms 



l-D s 



dtC{t) 



(22) 



For the further use together with the mode-coupling ap- 
proximations, the force-force correlator in eq. (22 ) should 



be rewritten in terms of the irreducible Smoluchowski op- 
erator 



(23) 



following [551 US] • After changing to the Laplace space 
according to C(z) = J °° dte~ zt C(t), 



C{z) 



(Ff 



1 



-F z 



(24) 



z-nr 

and using the standard operator identity for A = Ai+A 2 
(z - A)- 1 = (z- Ai)- 1 + (z- A)- 1 A 2 (z - At)- 1 , (25) 



(with A = fit, Ax = n tin ) for the resolvent in (24), one 
obtains the expression 



C(z) 



C'"(z) 



l + D s C**(z) 



(26) 



for the force-force correlator in terms of the irreducible 
one C l "(z) = (F z (z - & i ")- 1 F z ). Exploiting the rela- 
tion ( 26 ) for z = leads us to the desired result 



M = 



Ms 



1 + D S f^°dtC^(t) 



(27) 
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for the tracer mobility in terms of the irreducible tracer 
force autocorrelation function C iTI (t), which in the time 
domain is given by 



C'"(t) = (F:e< 



(28) 



Eq. 



(27) allows a simple interpretation if one intro- 

The friction 



duces the tracer friction coefficient ( 
coefficient is given by the sum 

C = Co + AC, (29) 

POO 

AC = / dtC'"{t), (30) 
Jo 

of the "bare" (single-particle) tracer friction coefficient 
due to the solvent, given by Co = M7 1 an( i the increment 
AC due to interactions with the bath particles. Because 
we are interested in dense dispersions where the fric- 
tion is highly increased beyond the solvent one, Eq. ( [29] ) 
provides a more secure route to approximations than 
Eq. (22 1. In Eq. (29) slow force fluctuations contribute 



to an increased friction, in Eq. (22) instantaneous and 



retarded velocity fluctuations need to cancel in order to 
yield a reduced mobility. The mode coupling approxima- 
tions to be performed later set up a self-consistent set of 
equations for slow fluctuations which is better tailored to 
Eq. ( 29 1 than to Eq. ( 22 ) . Recent simulations of a forced 



probe in a bath of noninteracting bath particles support 
to approximate C 1 " instead of C [27]. Even though the 
bath particles do not interact among themselves, the col- 
lisions with the probe particle induce correlations in the 
velocity but not (or to lesser extent) in the force fluc- 
tuations. The friction increment increases linearly with 
bath density, as expected from independent collisions of 
the non-interacting bath particles with the probe. This 
expected behavior, however, does not hold for the mobil- 
ity change, which varies more rapidly with bath density 

m- 

After these formally exact manipulations, approxima- 
tions are now required in order to evaluate the irreducible 
force correlation function. Low density approximations 
have been performed |28j , and we will use mode coupling 
approximations to address high packing fractions close to 
the colloidal glass transition. 



D. Transient tracer density fluctuations 

The important quantities in the mode-coupling ap- 
proach are the tracer and the bath densities 



p s (r) = 8(t-t.), 



N 



(31) 
(32) 



With the convention for the Fourier transform of a 
function X(r) to be 



(33) 



X(q) = / dre l * r X(r), 



implying 



X(r)= (2^/ dqP " q ' rX(q) 



for the back-transform, we have 



N 



pi = 



(34) 

(35) 
(36) 



for the tracer and the bath density modes. 



1. General properties 

The first question to clarify concerns the time- 
dependent correlator 



<P q e nt V q <) 



(37) 



of two tracer density modes. For which pairs of wavevec- 
tors q, q' is it nonzero ? 

For the case of an isolated system, translational invari- 
ance implies that after all particle positions have been 
shifted 



r -> r' 



(38) 



with r, 



a, i\- ->■ 



rj +a (i = 1, . . . , N), the average 
(|37]) should remain the same. Since the shift introduces 



the prefactor e 4 ( q+q • ) a in eq. (37) and the vector a can 
be arbitrary, the condition 

q = -q' (39) 

follows. 

Our driven system is translationally invariant as well, 
since we assume the external force F ex to be space and 
time independent and the Smoluchowski operator ( 14 ) 



does not change after the shift (38). Thus, the argumen- 



tation used for the isolated systems and thus the condi- 



tion (39) remains. 



We can thus introduce the usual notation 

K(t) = «e nt V q ), 



,(*) = 



NS, 



(p q e nt V q ) 



(40) 
(41) 



for the tracer and the bath density mode correlators, 



where S q is the bath static structure factor 



S i = N 



The Fourier back transform of 4>t[t): 
s (r)(t) = FT- 1 [^ q (t)]. 



(42) 



(43) 
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is the probability density of finding the tracer at the point 
r in space at time t with the initial condition that at time 
t = it was localized at the origin. The derivation is 
given in Ref. [55] > where an isolated system is considered. 
Since no special properties of the Smoluchowski operator 
for isolated systems are used in Ref. [29], the derivation 
is valid also for our case. 

The general condition for the Fourier back transform 



( 43 1 to be real is 



#-«(*) = #*(*)■ 



(44) 



We can easily see that the property ( 44 1 is indeed fulfilled 
by <£® if we use it's definition (40) and the fact that the 

operator e fit * is linear and real (since it contains deriva- 
tives with respect to phase space coordinates multiplied 
by real numbers). So, 



{e~ l q r = e° * e l q r = ) = <p° (i) (45) 



2. Zwanzig-Mom equations 

The Zwanzig-Mori projector operator formalism allows 
one to express the fluctuations of a given observable in 
terms of the corresponding memory kernel. After intro- 
ducing the projectors P A = A) (A*, Q A = 1 — P A for 
an observable A (we assume that (A* A) = 1 for simplic- 
ity), the Laplace transform C(z) = J °° dt e~ zt CA(t) of its 

equilibrium time correlation function Ca (t) — (A*e n>t A) 
can be written as [8] 



C A (z) 



z + lua- Ma{zY 



(46) 



where the frequency Ea and the memory function M A (z) 
are given by 



loa = -(A'ftfA), 
M A (z) = (A*tfQ A 



z - Qa^Qa 



(47) 

-Qa^A). (48) 



In the time domain, eq. ( 46 1 corresponds to 

dt C A (t) = -uiA C A {t) + [ dt'M A (t - t') C A (t'). (49) 
Jo 

For dissipative systems, like for the case of our system 
described by the Smoluchowski operator, a second pro- 
jection step is needed [25[j26]. To this end, one introduces 
the irreducible Smoluchowski operator 

fit'" = Qa (fi f - StfA)(A*StfA)- 1 (A*n J t) Q A . (50) 

and gets the representation 

M A (z) 



M A (z) = 



(51) 



for the memory function in terms of the irreducible mem- 
ory function M A (z) given by 



Mf(z) = (A'fit Qj 



1 



: Q A rtA), (52) 



which time evolution is generated by fit'": M A {t) 



(A*ti* Q A e nii,It Q A rtA). 
equation, one obtains in the time domain 



For the irreducible memory 



d t C A (t) = -u A C A (t) 



1 



W A Jo 



dt'M~ A (t-t')d' t C A {t'). 

(53) 

The procedure of expressing the equilibrium correla- 
tion functions in terms of memory kernels sketched above, 
was originally proposed for an isolated system evolving 
with the unperturbed operator fij. We apply it to the 
correlators evolving with the nonhermitian operator fit 
even though the mathematical conditions and justifica- 
tions are unknown at present. This procedure is based 
on the conjecture that the algebraic structure of the 
Zwanzig-Mori equations together with the mode coupling 
approximations, necessary in the latter steps to evaluate 
them, capture the mathematical bifurcation describing 
the derealization of the probe under strong force. At 
present this conjecture can only be tested by formulat- 
ing the theory and considering its results in comparisons 
with data and formal symmetry requirements. 

After these considerations, we are in a position to write 
down the memory equation for (f>i(t) = (p q *e nt *p q ) : 



a^*(i) = -w q 0'(t)- 

with 



dt'M^{t~t')d t ,^{t'), 
(54) 



w q = (55) 

M^'(t) = (p q *fitQ s e nti "*Q s fitp q ) (56) 

fit - = Q s (fit _ fit ) Wq 1 fit ) Q s (5?) 

and the projector 

Q s = l-p q )«. (58) 
Applying fit to p q yields 

fit Pq = D s {d 2 s + (F s + F ex ) ■ d.) e 4qr = (59) 
= D s {-q 2 + J q-(F s +F ex )) /9q , 

so for the frequency w q we obtain 

w q = D s {q 2 -iq-F ex ). (60) 

The fact that w q = — (p q | fit | p q ) turns out to be com- 
plex is obviously the consequence of the mentioned non- 
hermiticity of the operator fit = fij + D s F ex ■ d s with 

respect to the equilibrium- weighted scalar product Q. 
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We would like to discuss now the relationship between 



the irreducible operator introduced in eq. ( 23 1 for the 
force- force correlator (call it '"(i 7 ^ 2 )) and the one in- 
troduced in eq. ( 57 ) for the tracer density modes (call it 
^ G°q))- To this end, we set F ox = 0, go to the limit 
q — > in eq. (57) and employ relations (591, (60). Then, 



we readily see that 



into products of the two-point ones [T5] : 

= 5 k ,v W 4>L(t) NS p cf> p (t). 
Note also that as part of the approximation, the irre- 



lim n tirr (p*) = n tlr '(F s z ) (at F cx = 0) (61) ducible operator on the left hand side of eq. d67b was 



q->0 

holds. It is not surprising, since the (q — > 0, z — > 0) limit 
for the tracer density fluctuations is related to the tracer 
diffusion and in the absence of the external force, the 
well-known relation (3U1 GH] 



Di 



D, 



1 + lim, 



q->0,z- 



>0 (M q ' ,rr (q,z)/w q ) 



(62) 



for the long-time tracer diffusion coefficient implies 
the Einstein relation 



Di 



/' 



(at F P , 



0) 



(63) 



connecting with the long time tracer mobility p con- 
sidered in the last section. We see that the irreducible 
memory function plays the role of a generalized friction 
kernel. 



E. Mode-coupling approximations 

1. The memory function 

In order to obtain a self-consistent equation for from 
the memory equation (54), the irreducible memory func- 



tion ( 56 ) is treated using the standard approximation 



steps of the mode-coupling theory |16j . To this end, first, 
the projectors onto the space spanned by the tracer-bath 
density products 

Pi= £ ^ P ). 9 (k,p,k',p')(p k ^ (64) 

k,p,k',p' 

are introduced, where the normalization matrix g has to 
obey the condition 

£ (Pk'iV I PwPp) 5(n, m, k', p') = <S n , k <5 m:P (65) 
k', P ' 

which requires an arbitrary vector | p k p P ) from the space 
of the tracer and bath density products to be left invari- 
ant upon application of P 2 S . 

As the first approximation, the "fluctuating forces" 
Q s W /Oq in eq. (56) are replaced by the projected ones: 



M* 



P° 



(66) 



In the next approximation step, the fou r-po int correlators 



(Pk Pp efl 1 Pk'Pp') appearing in eq. (66) are factorized 



replaced by the normal one on the right hand side. 

It is easy to see that for t = 0, the factorization ap- 
proximation (67) becomes exact and allows us to calcu- 



late the normalization matrix g from the condition ( 65 1 , 
which then reads 



£ £k',k S p > >p NS p g(n,m,k',p') = S n ^S m>p , 

k',p' 



leading to the result 

,g(n,m,k, p) 



NS V 



(68) 



(69) 



The only missing parts appearing afte r th e application 
of the projectors P% in the expression (66) are now the 
static averages of the form 



(70) 



and 



(pgrtQ'pfr,,). (7i) 

First, we notice that since is not self-adjoint, the 
are not complex conjugated of each 
be the case for an isolated system. 



averages (70), (71 



other, as it wouli 



Using eqs. (58), (59), for the term (70) we obtain 



D 



UpuPpQ s ^pD = {p 



q/ - \f q-k P*p(-1 + i q • ( F * + F e 

(Pq-kP P X-9 2 + *q-(F s +F e:c )). 



■M 



(72) 



The nontrivial terms in the above expression are either 
of the form of the tracer-bath static structure factor 



(73) 



or the term of the form (F s p k p p ). The latter can be re- 
duced to the tracer-bath static structure factor by means 
of partial integration 

(FsPLp P ) = ^ J dT(d s e- v )plp p = 

= -~J dVe- v d s {pl Pp ) = -ik6 k+pfi S° p . (74) 

Using this relation, one readily obtains the result 

(pi*p* p Q s & Pq > = <5 q , k+p D s S s p (q ■ p) (75) 
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To calculate the remaining average (71), the fastest 
way is to we use the adjoint of 



<j%rtQ'rij> 9 ) = ($r*p'Z)Q'Ap v ) 



(76) 



and its action on the conjugated tracer density mode (see 

e q . my. 



n ad 3p s * = D s (-q 2 -iq-(F s -F ex )-F ex -F s )p s *. (77) 



The calculations go in the same way as for the term ( 70 ) 



and we omit the details except for the fact that compared 
to the expression (59), an additional term F ex ■ F s is 



present in the brackets in eq. (77), so that the final result 
also contains F„ x : 



^ +p D s S s p (q-p-iF ex -p). (78) 



After collecting the terms together, using relations 



(66), (67), (69), (75) and (78), the final mode-coupling 



expression for the memory function reads 



M«''"(t) 



k+p=q 



Pis; 2 



q p (q p - i F ex -p) (f>l(t)(f> p (t). 

(79) 



2. The bulk dynamics 

So far, nothing was said about the fluctuations of the 
bulk density modes (/>q(t). The simplest reasonable as- 
sumption for these is, to neglect the effect of the external 
force. In the thermodynamic limit, which is considered in 
this article, this assumption is justified, since the effect 
of the tracer on the bath will be to perturb its neigh- 
bourhood only locally. This effect will be included in 
our theory and can be described by the tracer and bath 
density mode products p q /0q, which are the Fourier space 
counterparts of the relative bath-tracer density p(r — r s ). 

So, we assume the bath dynamics to be unaffected by 
the external force. The memory equation for the bulk 
dynamics reads 

r 9 cW q (i)+0 q (t)+ / dt'M™(t-t')d t 'cf>a(t') = (80) 



with 



Tg = S q /(D q 2 ) 



(81) 



Results from these well studied MCT equations jTS], [TB] 
will be used in the following whenever properties of the 
unperturbed bath particles are required. As most impor- 
tant result let us recall already here that MCT predicts a 
glass transition of colloidal dispersions at high concentra- 
tions when density fluctuations do not relax completely. 



3. The force-force correlator 

In order to obtain the mode-coupling approximation 
for the irreducible force-force correlation function ( 28 1 , 
which enters the expression ( 29 ) for the tracer friction 



The standard equilibrium MCT expression for the ir- 
reducible bulk memory function is 

M q "(i) = ^4 E nS q S k S p (q • (kc fc + pc p )) 2 MWpit) 

k+p=q 

(82) 

where n = N/V denotes the number density of the bath 
particles and c q is the Ornstein-Zernike direct correlation 
function: 



coefficient, we use the usual MCT strategy and substitute 
the force F GX by the one projected to the tracer-bath pair 
density modes: 



c""(t) = (F/ e ° tirrt F/)«(F;p 2 s 



tp » F , 



(84) 



where P| is given by eq. (64) 



Eq. ( 84 ) is the simplest possible approximation of the 



mode-coupling type for the tracer force autocorrelator, 
since a product of at least one bath and one tracer den- 
sity mode is needed. This is due to the fact that the equi- 
librium average with the tracer force (F* . . .) appearing 
in eq. ( 84 ) is zero both for a single tracer and for a single 



bath density mode. 

The calculation is completely analogous to that of the 
memory function in the previous paragraph. One uses 
the factorization approximation (67) as well as the rela- 



tions (69), (731 and (74) to obtain 



C'"(t) 



k 



1 

NSk 



k, St, 



(85) 



F. Reality of the observable averages 

In order to check, whether our MCT approximations 
preserve the reality of observable quantities, we choose 
the friction coefficient as the typical example. The MCT 
expression ( 85 1 for the irreducible force- force correlator 
C iTT (t) enters the eq. (30) for the friction coefficient in- 



cremen t un der the time integral. 

Eq. (85) contains the sum (over all k £ K 3 ) over 0£ 
multiplied with real and rotationally invariant (in the k- 
space) factors. The similar structure arises also if one 
makes mode-coupling approximations for other (tracer- 
related) correlators, since the dynamic part is given by 
the factorized four-point tracer-bath correlators and the 
response quantity-specific part comes in via the different 
static k-dependent "vertices" . 

In order for the friction coefficient to be real, it suf- 
fices to show that the MCT approximation for the tracer 
density correlator fulfills the condition ( 44 ) since then 
the imaginary parts in the k-sum cancel. 



S q = 1/(1 - nC q ). 



In the mode-coupling equation (|54|) for q , the memory 



(83) function ( 79 1 couples the correlatofior the wave vector q 



to the correlators for all the other wave vectors, one has 



to consider eq. ( 54 ) as a system of coupled equations for 
the set of all q. 



A rigorous proof that (44 1 will hold for t > 0, pro- 
vided that the initial conditions (l86|) hold, however, was 



not given here. Such a proof would require the explicit 
construction of the solution. The numerical results for 



From a purely mathematical standpoint, it can have 
different solutions depending on the initial values </> q (t = hard sphere glasses in Sect. IIIB violate (44) for strong 
0) which in general might not fulfill the condition (44). 
But in our physical problem 



K(t = 0) = 1 



(86) 



holds so that (44) is fulfilled for the initial values. 



We show now that the assumption that ( 44 ) is valid for 



t > does not contradict the system of equations ( 54 ) 
For this purpose we look at the equation for — q: 



where 
and 



-c_ q </>> (t) - / dt>Mt' l ;*(t-t') <wi q (0, 

Jo 

(87) 



D s (q 2 + iq-F ex ) 



E 



forces in a narrow angle of directions around the one per- 
pendicular to the external force. This will be discussed 
in more detail in the later section and could indicate a 
break-down of the theory for very large forces. Yet, we 
continue on the assumption that our qualitative results 
are not affected. 



III. RESULTS FOR THE HARD SPHERE 
SYSTEM 

Here we apply the microscopic formalism developed 
in the last Section, to the colloidal hard sphere system. 
The tracer radius a s can be different from that of the 
bath particles a (see Figure [I]). The control parameters 
are then the ratio a — a s /a the radii of the tracer and 
the bath particles and the volume fraction of the bath 
particles ip 



4 3 

I ira n. 



q-p' (q-p'+iF^-p') p ,(t) 



Assumption (44) yields: 

M^'(i) = (M q -T(*)- 



(89) 



(90) 



A. Low density limit 

First, we would like to calculate the tracer friction co- 
efficient increment A£ in the limit cj> —> of vanishing 
volume fraction of the bath particles. We use the rela- 



To prove this, we notice that in the expression ([79} for tion (30j) for AC in terms of the force correlator C i,r (t) 



the memory function every term for a certain pair of wave 
vectors (k, p) is complex conjugated with the term in 
expression ( 89 ) , corresponding to the pair of wave vectors 



together with the MCT approximation ( 85 1 for C"' (t) 
and obtain 



(k', p') with k' = -k, p' = -p: 

q • (-p) (q ■ (-p) + iFe* • (-p)) 01 k (t) 0- P (i) 

~rs2 



AC 



1 



(2tt) e 



<3fe 



dt&(t)Mt)- (93) 



Note that the k-sum in ( 85 ) was changed to the integral 
over the k-space: J2k ~~ ^ V/i^n) 3 J dk 



S s / q-p (q-p + iF B!B -p)CWp(t) 



The static structure factors Si, Sk entering eq. (93) 



/ 1 



Sf q'P(q-p-»F a -p)^(t)^(t) 



\NS P p 



(91) 



can be easily calculated to the leading order in 
start with the tracer-bath structure factor: 



So, given that (44) holds, (90) holds also. On the other 
side, if we use (44) in eq. (87), we get: 



S' h = 



N 



ik-(r 



rs) ) fa Y(e lk ' (ri " rs) ). 



We 



(94) 



(92) 

and this equation is the complex conjugated of eq. (54) 
for </> q due to relations (88) and (90). 



The approximation made here consists in considering the 
system as a conglomerate of independent two-particle 
clusters. The two-particle structure factor (e lk •( ri ~ rs )) 
can be reduced to an elementary integral over M 3 and 
yields a Bessel function: 



These considerations snow thaxme condition ( 44 ) is 



consistent with the mode-coupling equations ( 54 ) , so that 



Si. 



-A-Knd A 



we can state that the condition for the reality of the 
Fourier back transform at least can be imposed on the set 
of mode-coupling equations for the tracer density mode 
correlators. 



suxx 



COS X 



(95) 



where x = kd and d = a. 



a is the sum of the tracer and 
bath radii. Eq. (95) shows that Sj* scales as Sj* = 0(<j>) 
with d>. 
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As for the bath structure factor, it can be written as 



S k 



1 



(96) 



for the case that the tracer is identical with the bath par- 
ticles. This follows directly from the definition of Sk (see 
eq. 42) since all the particles in the system are equiv- 
alent. So, the bath structure factor scales with <j> as 
Sk = 1 + 0{4>). Since the structure factors S% and Sk 
enter as a product into the relation ( 93 ) , to the leading 
order in (j), the expression ( 95 1 for Si 
0-the order approximation 



together with the 



(97) 



can be used. 

For the time-dependent density fluctuations </>£(£), 
<f>k(t), the 0-the order terms in <p are non-zero and 
can be obtained by neglecting the memory integrals in 
eqs. ([SlLpO} and setting Sk = 1 in eq. ((811): 



ik(*) 



a D s (-fc^+ik-F cx )t 



-D k t 



(98) 
(99) 



So, to the leading order in <f> the higher-order corrections 



to </> k (i), 0k (t) can also be neglected in eq. (93). 



Note that according to the translation theorem of the 
Fourier transform theory, </> k (t) in eq. (98) is the Fourier 



transform of the Gaussian, whose maximum position 
shifts by the distance D s F ex t from the origin with the 
time t. It is not surprising, since we made the approx- 
imation of completely neglecting the effect of the bath 
particles, i. e. considering a single tracer particle in a sol- 
vent, scdimenting under the external force F ex so that 
it's mean position at time t is given by D s F ex t. 



With the approximations (98), (99), the time integral 



in eq. ( 93 ) can be performed with the result 



1 



(D + D s )k 2 -iD s k-F ex 
(D Q + D s )k 2 +iD s k-F ex 
(D + D s yk* + D*(k-F ex )* 



(100) 



We notice that the imaginary part in the above expres- 



sion does not contribute to the k-integral in ( 93 ) , since 
it is antisymmetric in k, so that the result for A£ is real. 

We choose now the k z -&xis in the direction of the ex- 
ternal force F ex and change to spherical coordinates. 



Eqs. (93), (100 1 then yield AC in the form 



with the dimensionless friction increment function 



(101) 



18 



dx 



SIM 



— COS X 



.arctan^ 
(102) 
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FIG. 2: Low density results for the dimensionless friction 
coefficient increment as function of the dimensionless external 
force parameter j3 = F cx a/(fcsT). Continuous line: exact 
result from Ref. [25] ■ Dashed line: our MCT result. 



and the dimensionless external force parameter 



P = 



D s {a 



)F K: 



D +D s 



(103) 



After employing the Stokes-Einstein relation ^ for Do, 
D s and reintroducing the factor fcgT the parameter /3 
can be expressed as 



F ex ci 

~k~wr 



(104) 



and has thus the physical meaning of the work done by 
the external force over the distance of bath particle radius 
in units of fc^T. 



We compare now our result ( 101 ) for A£ with the exact 
low-density result from Ref. [28 , where the two-particle 
Smoluchowski equation was solved for arbitrary exter- 
nal force values. The differences appear in the numerical 
prefactor (1/2 in the exact theory instead of 1/3 in our 
calculation) and the dimensionless friction coefficient in- 
crement functi on (s ee Fig. [5]). Despite these differences, 
the expression (104) for the dimensionless external force 



parameter and the scaling of the friction coefficient in- 
crement with the system parameters <p and a, agree with 
the exact theory. 

From Fig. [2] we see that for values of (3 ^ 2, where the 
initial decay from the linear response plateau occurs, the 
MCT result agrees well with the exact low-density cal- 
culation, whereas for values of /3 ~ 10 the decrease of 
the tracer friction coefficient with F ex is strongly overes- 
timated by the present version of MCT. For f3 — > 00 the 
exact theory predicts a second plateau value of 1/2 for 
f{f3), which is zero in our MCT calculation. 

Note that the approximation made to obtain the ex- 
pression ( 98 ) was to completely neglect the memory inte- 
gral, and the latter also contains the external force. Such 
an approximation might not be valid for large F ex . The 
large /3-limit turns out to be singular in the low density 
theory [55] and this could also hold for our mode-coupling 
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equations. This issue surely deserves a more thorough 
analysis including the numerical solution of the full time- 
dependent equation for the tracer density correlator. 

The elementary consideration here is just aimed to 
show that even the simplest reasonable approximation 
(981 for the tracer correlator, in combination with the 



mode-coupling expression ( 85 ) for the force-force corre- 



lator, is already capable to explain the thinning effect. 
This is due to increasingly strong time oscillations of the 
integrand in eq. ( 93 ) , which cut off larger and larger por- 



tions of the integrand with increasing F e , 



B. Long time limit of the tracer density correlator 

After considering the low-density dynamics in the last 
section, we turn now to the high-density "statics". I.e. 
we will assume that we are above the critical point where 
the colloidal bath forms a glass (ip > ip c ). The first ques- 
tion then concerns the existence and force dependence of 
the long time limits 



= lim ^(t) 
L = lim q (t) 



(105) 
(106) 



of the tracer and the bath density correlators. We restrict 
ourself here to the case of the equal size of the tracer 
and the bath particle so that a = 1 holds and only the 
bath volume fraction tp remains as the control parameter. 
The existence of a finite probe 'non-ergodicity' parameter 
signals that probe density fluctuations remain frozen-in, 
and that the probe can not explore the whole volume V. 
A /' > signals that the probe is localized. 

For t — > oo, the full time-dependent MCT equations 
(54), (80) reduce to the "static" (i. e. time-independent) 



equations 



/q 



q - M^"(r, /), 
Mq"(/), 



(107) 
(108) 



(the derivation of the eq. (107) for the tracer straight 



forwardly follows the steps of the standard MCT |16j). 
Equations (107), (108) are much easier to treat than the 



full time-dependent equations because no memory inte- 
grals over many time decades have to be evaluated. Still, 
they contain a wave vector integration and hence repre- 
sent a coupled system of equations for all q. 

In order to be solved numerically, the long time limit 



equations ( 107 ), ( 108 1 have to be rewritten as fixed point 
equations: 



/q 

/a 



M" 



(109) 
(110) 



Starting with the initial values / q — 1, /* = 1 for all 
q, the right-hand side of eqs. (110), (109) is evaluated 
and then used as input for the next iteration, until the 
desired precision is reached . 

The spatial symmetry properties of the tracer and bath 
correlators dictate the choice of the discretization scheme 
for and / q in the g-space. Since the bath correlator 
is not affected by the external force, as was discussed in 



Section HE 2 f q retains its spherical symmetry so that 
it can be taken as input for the equation (110) from the 
standard isotropic MCT calculation. An equidistant q- 
grid is chosen here with 100 g-points for qa G (0, 20) 
and the so-called "Bengtzelius trick" (see [32 for details) 
allows to reduce the computational complexity during the 
calculation of the memory integral. 

As for the external force introduces a preferred di- 
rection in space and thus breaks the spherical symmetry. 
This means that in contrast to the isotropic calculation in 
the absence of F ex , one explicitly has to resolve the angu- 
lar structure of / q . One symmetry still remains, namely 
the rotational symmetry around F ex . Thus, it is enough 
to specify the magnitude q and the angle 8 between the 
vectors q and F ex in order to determine / q uniquely. 

We thus introduce spherical coordinates with the z- 
axis pointing in F cx -direction and assume the wave vector 
q to lie in the x-z plane. The property /£_ = /i* (which 



IIF) 



we assume to be valid, see discussion in Section 
allows us to confine the angle 8 to the interval (0 
Note however that despite this quasi-2d nature of the 
problem, the integration in the memory kernel includes 
all the wave vectors and cannot be reduced to a two- 
dimensional one. So, the polar angle <p also has to be 
included there, too. 

Since the angles 8 enter in eqs. (110) only via scalar 
products of vectors q, p and F cx with each other (see 

i. e. only as cos 6 



( 79|) for w q and 



relations (60), 

or sin 8, it appears reasonable to change the variable from 
8 to y — cos 8 and employ the relation J™ dd sin 8 f{8) — 

J dy f{y) in calculating the wave vector integral entering 



Since the / q -s enter the tracer equations ( 110 ) as input, 
the same equidistant g-grid is chosen for / q as for / q . The 
y- and </>-grids are also chosen equidistant with 30 points 
for y G (0, 1) and 20 points for <f> G (0, 2tt). 

Finally, as the equilibrium structural input, the 
Percus-Yevick structure factor [33l [34] is used. 



The results of our numerical calculations for ff, are 
n 

shown in Fig. k3l The value of the control parameter 
tp was choose to be 0.52, i. e. closely above the glass 
transition point cj) c w 0.516 of quiescent MCT. The real 
and the imaginary parts of / q are plotted as function 
of q and the different curves for the same force value 
correspond to the different values of angle 8. F, 
in the dimcnsionless units of ksT '/a (see eq 



, iX is given 
The 



104) 



magnitude of the wave vector q is given in units of 1/a. 

For F ex — 0, no angle dependence is observed and 
the standard isotropic MCT result is recovered. As F ex 
increases from and 30, the splitting of / q -s with 8 be- 
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FIG. 3: Long time limit of the tracer density correlator (upper panel shows the real and lower panel the imaginary parts) 
as function of the wave vector magnitude q for different angles 9 (see text) and different F ex . The bath volume fraction is 
ip — 0.52. The two /f at F ex = in the upper panel show that the used angle-discretization is fine enough to recover the known 
isotropic result. In the lower panel the collective f q is included as reference. The rapid variations in all f q below qa < 0.5 are 
discretization artifacts which do not affect the results for larger q. 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 




FIG. 4: 2d contour plots of / s (r) (equi-probability contours are labeled by their values of f s ) for if — 0.52 and F ex = 
0, 10, 20, 30. The external force points in the positive horizontal direction and r is given in units of the tracer radius a s — a. 
Data from Ref. 14 , where a smaller region was shown. 
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FIG. 5: Critical force vs. bath volume fraction. Data from 
Ref. |14j . The dashed line is a fit with eq. ( 131 1 for the 
schematic model, where the relationships f ~ f + Cy^e, 
e = (4> — 4>c)/4>c where assumed. 



comes more and more pronounced, revealing the spatial 
anisotropy, which increases with the external force. 

For F ex > 40, /q appears to drop to zero for all q. This 
is an indication of the bifurcation transition in the equa- 
tions (107): there is a certain critical force FiL, above 



which only the (always existing, trivial) zero solution of 
( 107 1 becomes stable, whereas below F£ x the numerically 



found non-zero solution is stable. 

Unfortunately, the numerics becomes unstable for 30 < 
F cx < 40, so that we cannot approach the region of the 
critical force value arbitrarily closely. But the observed 
trends in the F^-dependence of the curves suggest that 
the transition is continuous ("of type A" in the MCT 
classification). I. e. for every q, /* goes down to zero con- 
tinuously with increasing F ex . At the bifurcation point, 
/q becomes zero so that the two solution branches of 



eq. (107) coalesce. 

The physical meaning of the bifurcation transition is 
that the tracer depins from the glassy matrix, it gets 
"pulled free" out of the cage by the external force. 

If we perform the Fourier back-transform of our q- 
dependent data (see eqs. (33), 
on the Fourier transform): 



( 34 ) for our convention 



r(r) 



1 



(2tt) j 



(111) 



we obtain the probability density / s (r) of finding the 
tracer at space point r for t — ¥ oo after it was localized 
at the origin at t = 0, as was discussed in Section II D 1 



The results for f s (r) are shown in Fig. [4 as 2d-contour 
plots. They can be intuitively interpreted as the shape 
of the cage around the tracer. In the absence of the 
external force, there is no distinguished direction and the 
cage is isotropic. If the external force if applied, we see 
that the cage gets deformed and elongated in the force 
direction. The maximum of the distribution (the cage 



"center") moves in the direction of the force. For F ex = 
10, the maximum shifts from the origin by approximately 
0.15 a in the force direction. Also a slight anisotropy 
around the maximum position can be seen. 

For F ex = 20, 30, the maximum position lies at almost 
the same value of approximately 0.25 a in the force direc- 
tion from the origin. The anisotropy of the cage shape 
increases quite dramatically, so that the cage acquires 
the shape of an "egg" with its thick side showing in the 
force direction and slightly dented. One possible explana- 
tion for this subtle feature could be due to the increased 
density of the bath particles immediately in front of the 
tracer in the direction of the force, as compared to the 
other directions. We expect this on physical grounds, 
since the bath particles should get compressed by the 
tracer moving in the F ox -direction; yet calculations of 
the bath particle densities have not been performed. 

Note that f s (r) stays positive practically in the whole 
region (r/a) < 1 of r-points considered. Whereas the 
reality of f s (r) can at least be imposed easily, as it was 
actually done in our numerical calculation here, the posi- 
tivity of the Fourier back-transform of the given function 
is a non-trivial matter. The fact that our MCT approx- 
imations preserve the positivity over a considerable spa- 
tial range should be appreciated. 

The quantity / s (r) has a very clear physical meaning 
and the results for it, discussed above appear reasonable 
and intuitive. They were obtained by numerical calcula- 
tions from the numerical data for which themselves 
are not so easy to interpret at first glance. It is still pos- 
sible, however, to relate the observed features in the be- 
haviour of f s (r) (Fig.[4| to the features of the /q-curves 
(Fig- §■ 

First, we see that for F ox = 10, the real parts of /q(q) 
are monotonously decaying and the imaginary parts have 
one maximum for all values of 9. Furthermore, the 
shape of the curves resembles the one obtained from 
the 9- independent real function /q(F 0X = 0) by multi- 
plying by the prefactor exp(iq • Sr) = cos(q Sr cos 9) + 
i sin(q Sr cos 9): 



/n(F cx ) « / q (F ex = 0) exp(iq • 8t(F„)), 



(112) 



where we assume 5r||F cx . Due to the translation theo- 
rem of the Fourier transform theory, this prefactor corre- 
sponds to the shifting of the / s (r)(F cx = 0)-distribution 
by vector Sr ~ F cx . This is in accordance with our ob- 
servations for the behaviour of / s (r). 

For higher values of F cx , for some values of 9 the 
real parts of the curves /q(q) become non-monotonous, 
whereas the imaginary parts exhibit a maximum and a 
minimum. This complicated behaviour can be related to 
the subtle features in the shape of the / s (r)-distribution 
discussed above (anisotropy with respect to the maxi- 
mum, "dent"). 

Finally, we discuss some peculiar features in the be- 
haviour of the imaginary part of For the direction 
perpendicular to F cx , the imaginary part of /* should 
vanish. This is the consequence of the relation /£. = /** 
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and the rotational symmetry around F ex , which implies 
/q* = for q X ■ F ox = 0, so that 



»{/q}-0 for ql_F e 



(113) 



should apply. 

For F ex = 10, we indeed observe this behaviour (which 
is also in agreement with our approximation ( 112| )), 
whereas for higher force values, a strong dip is observed 
for 3{/q}(g) at qa ps 3, for the 0-values closest to tt/2, 
whereas at q = 0, 3{/q}(<?) seems to diverge for all 9. 
Since 9{/q} is antisymmetric in q, this behaviour sig- 
nalizes the discontinuity of the function 3{/q}(q) at the 
plane q_LF cx . 

It is easy to see that the origin of this discontinuity lies 
in the presence of the term w q , which enters the long time 
limit eq. (110) under the denominator on the right hand 
side. The term cj q = D s (q 2 — iq • F ex ) itself obviously 
has the property 3{w q } = for q±.F ex , so the same 
should be valid for its inverse. The inverse however, is 
divergent at the point q = 0, so that S{l/w q } = applies 
everywhere on the plane qUF ex except for the singular 
point q = 0. 

The discussed behaviour of the term l/w q is relevant 
for since if one approaches the bifurcation transition, 
the function /q(q) becomes confined to a narrow region 
around the origin in the q-space and the same applies 
also for the memory function Mq' irr . So, as a very rough 
first guess, i\^fq ,lri: can be neglected in the denominator 
of eq. (110), as a "small" quantity. As a next approxi- 



mation, one can use an q-independent ansatz for Afq ,llr 
in (110). In this way it is already possible to obtain the 



structure of the q-dependence very similar to that ob- 



served for f* in our numerical solution of eq. (|110| at 



high force values. 

We believe that the qualitative discussion of the prop- 
erties of given here, already contributes to a deeper 
understanding of our numerical results. The origin of the 
unphysical results for /* at larger forces remains to be 
clarified. It should be possible to extract also quantita- 
tive results for the limiting cases F cx — > and F ex — » F c c x 
by either performing expansions in F ex (for F cx — > 0) in 
eq. (110) or an appropriate asymptotic expansion assum- 
ing the smallness of (for F cx F£ K ). This will be a 



subject of future work. For the case of schematic models, 
the latter has recently been achieved [22] . 

Up to now, only the results for (p — 0.52 were dis- 
cussed. For higher volume fractions, the results for 
look similar, but the critical force increases with ip. The 
results for F c c x vs. ip are shown on Fig. [5] which can be 
considered as a phase diagram for the tracer states (im- 
mobile below the F c x (<p)-line and mobile above). At all 
densities, the numerics became unstable upon approach- 
ing F c c x . Error bars in Fig. [5] mark a region of confidence 
in the convergence of different numerical solution tech- 
niques, to estimate the critical forces. For the values of 
F cx outside the range marked by the error bars, the al- 
gorithm converged and produced a localized state of the 
tracer for F cx below the lower limit and a delocalized 
state for _F 0X above the upper limit. 



IV. SCHEMATIC MODELS 

A. Construction of the models 

Simplified schematic models for tracer density fluctu- 
ations can be constructed by considering only two wave 
vectors k, — k with k 1 1 F ex and neglecting the contribu- 
tion of all the other wave vectors to the memory function. 
Due to the property 01 k = </> k * , the system of equations 
for </> k , 4>-u> obtained in this way, reduces to the following 
single equation for ^> k : 



d t &(t) = -D s {k 2 -ikF ex )<j> s k (t) 



(114) 



- / dt'v^^t-t' f ^ k (i-t')d t ^(i'), 

Jo 



where Wk = 4S"^ /iVS^fc. We see that th e for ce depen- 
dence of the memory function drops in eq. (114). Further 



simplification can be achieved if we introduce a scale for 
lengths so that k = 1, rescale time so that D s — 1 and 
model the behaviour of the bath correlator (j>k=i = 4> 
with the standard F12-model [15]. For simplicity we set 
its short time coefficient equal to the tracer's. 
We finally arrive at 



d t 4> s {t) = -u>(p s {t)- / dt' m s {t-t')d t i4> s {t'), 



dt<Kt) = -<P{t)~ / dt'mit-t^dpW), 
Jo 



(115) 
(116) 



with the complex frequency 



to = 1 — i F e 



(117) 
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and the memory functions 



m s (t) = v s cp s *{t)^{t) (118) 
m(t) = Vi<t>(t) +v 2 cj) 2 {t). (119) 
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FIG. 6: Real and imaginary parts of the long time limit of 
the tracer density correlator as functions of the external force 
(v 2 = 2.0, v a = 4.0). 



Considering the real and imaginary of eq. (121 1 sepa- 
rately, we obtain the system of equations 



*i(A*,/ 2 *) = o 
j ? 2 (/r,/i) = o 



for /f , /|, where 



F x = ft + Fevfi - v.f (ft - ft 
F 2 = fS-F aB f(+v.ffi. 



Eq. ( 123 ) yields the linear relationship 



between /J and /* with 



(122) 
(123) 



fi ) (124) 
(125) 



(126) 



This schematic model will be called the F ex -Sjdgren 
model. It extends the original Sjogren model (proposed 
in Ref. [35]) for the tagged particle dynamics in a host 
fluid for the case of external driving with the force F ex . 
From the described rescalings, it is obvious that in the 
schematic model F cx is measured in thermal energy di- 
vided by the length used to rescale the wavevector in 
Eq. JITih. 



A modification of the model ( 115 ) is of interest, namely 
the one where cf>(t) = 1 is set: 



d t cf) s (t) = -uj(f> s (t)- dt'v s (t> s *(t-t r )dt><j) s (t') . (120) 



1 



f 



(127) 



which can be inserted into eq. ( 122 1 to obtain a quadratic 
equation for Its first trivial solution is /f = 0. The 
second solution is: 



fl = 



-a F K . 



f-1 



/(1 + a 2 ) 



(128) 



In the absence of the external force we have a — 0, so 
/| = 0. I. e. the tracer correlator is real and the result 
of the original Sjogren model is recovered: 



ft 



1 



1 

Vs f 



(129) 



We call this model the F ex -Fl model. It extends the There is a bifurcation at v c s = l/f: for v s < v%, the long 
well-known Fl-model [T5] (corresponding to the case 
Fex — 0) which was introduced in connection with the 
mode-coupling theory for the Lorentz model. The lat- 
ter describes a tagged particle in an array of immo- 
bile randomly distributed scatterers. As we shall see in 
Sec. IV G 1 the z — value for the Laplace transform of 



time limit of <j) s {t) is zero (since it cannot be negative) 
and for v s > it has a non-zero value which is given by 
eq. (129). This consideration is valid for / > 0, i. e., if 
the bath is in a glassy state. If / = 0, then /* = always 
holds. 

Now, let us assume that for F ex = 0, the tracer is in 
an arrested state: ff > 0. Can this state be "molten" 
by the external force field ? We look at the behaviour of 
fi as a function of the external force for fixed values of 
/ and v s . Eqs. (p7|,(p8| yield: 



B. The phase diagrams 

In this section we consider the long time limit / s = 
f 1+1/2 of the correlator <fi s (t). For our schematic models, 
the condition (107) reduces to 



l-f s 



m s (f s J), 



(121) 



f s _ (l+Vsf)((Vsf) 2 -l-Fl) 

In Fig. [6j /* and /J are plotted as function of F ox for 
fixed parameter values v 2 — 2.0, v s = 4, . We see that 
our schematic model exhibits the critical force at which 
the tracer becomes delocalized, i. e. where the depinning 
transition occurs. 
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mobile / localized 



1 /f 




FIG. 7: Phase diagram sketches of the schematic model for 
fixed f > f c (upper panel) and fixed v s (lower panel). 



From eq. (130), we obtain 



(131) 



for the critical force F£ x , at which real and imaginary 
part of f s become zero together. Eq. ( |131 ) can be visu- 
alized in form of phase diagrams, (schematically) shown 
in Fig. [7] Note that in the fixed v s diagram the horizontal 
axis corresponds to the variable e(f), which measures the 
distance from the glass transition point of the bath (see 
Sec. IV D for the definition) so that s = corresponds 
to / = f c . This allows a (qualitative) comparison to 
the microscopic phase diagram in Fig. [5j where we fit- 
ted the p oints from the g-dependent calculations using 
Eq. (131) and assuming the relationships f = f c + C\fe 
and e = (<f> - (f> c )/(f) c 



The latter lead to the expression 
1 (with a = v s f c and 



(0 + (b/y/0~ c )y/(j)- 4> c ) 2 

b = v s C) which was used for the fit. The values of the fit 
parameters were a — 4.36 and b — 327.61. 

The analytical results for the long time limit derived 
above for the -F^- Sjogren model obviously apply for the 
F ex -Fl model if one sets / = 1 in the corresponding 
expressions. 



C. Bifurcation analysis 

The purely algebraic consideration in the last Section 
can be completed by the analysis of the bifurcation sce- 



At the critical force, the two solutions f s = and the 
one given by eqs. (128), (126) coalesce. In order to clar- 



ify the geometry of the problem, we look at the solution 
space of eq. (121), namely the (/-f, /|)-plane. Since the 
complex conjugation is involved, being a nonanalytic op- 
eration, no use can be made of the complex analysis |38j . 
So, instead of the complex eq. (121) we have to ana- 



lyze the equivalent system of real equations ( 122 ), ( 123 1 



Each of them defines a curve in the solution plane. The 
curves are of the order not higher then quadratic. Gener- 
ically, they intersect in two points. The bifurcation oc- 
curs, when the two curves are just tangent to each other, 
so that there is only one intersection point identical with 
the osculation point. 

The normals to the curves defined by eqs. ( 122 ) 
and ( 123 ) are given by the gradient vectors 
(dF 1 /8j[7dF 1 /dft) and (dF 2 /dft, dF 2 /dfi), re- 
spectively. They have to be parallel at the bifurcation 
point (since the curves have to be tangent, as discussed 
above) and this is equivalent to the requirement 



det J = 



(132) 



for the matrix J with the elements Jik = dFi /<3/| , where 



, fc G {1,2}. From eqs. (124 1, (125) we obtain 



with 



I = 



J = 1-1 



v s f -F s 



F, 



-V s f 



(133) 



(134) 



and the unity matrix I. So, condition ( 132 ) is equivalent 
to the condition det I = 1, which immediately yields the 



expression ( 131 1 for the critical force. 



Besides from recovering the result for the critical force 
from the last section, the considerations here enable us 
to learn more about the bifurcation. Let us look at the 
eigenvalues A of the matrix I. They are solutions of the 
characteristic equation det(I — A 1) = 0, which yields 



A 2 = vlf 



F 



(135) 



The eigenvalue A = 1 corresponds exactly to F ex = F£ x . 
This eigenvalue is not degenerated. This means, the bi- 
furcation is of the codimension one 1 36 . 



D. Effect of the external force on the correlators 

We want to consider now the full time dependence of 
the correlators and look, how it is influenced by the ex- 
ternal force. The equations (115), ( |120[ ) are solved nu- 
merically using the algorithm described in [37] applied 



to real and imaginary part individually. The results are 
summarized in Fig.[8]for typical values of the parameters. 

The behaviour of the schematic models of the standard 
MCT well known from Refs. QH QH |35] will be taken as 
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reference where appropriate in the following. For the 
parameters of the F12 model the relation holds 



(- 



1) 



(136) 



for Vi, vl lying on the bifurcation line, separating the 
liquid phase from the glass phase. We introduce also 
the parameter e which measures the distance from the 
glass transition line and is related to vi = v% + <5i»2 and 
Vi = + dvi by 



Svi f c + 6v 2 r 



(137) 



with f c = 1 - 



In the following, the parameters will 



be chosen such that v-i = v%, so that <5t>2 = 0. Thus, spec- 
ifying the values of vi and e completely determines the 
parameters of the F12-model according to the eqs. ( 1371, 
p6|. 



First we consider the case, where the tracer is in the 



arrested state for F c . 



0. For the F ex -Fl model this 



means v s > 1, and for the i^^-Sjogren model e > 
and v s > 1//. The corresponding plots are shown on 
panels (a,b) (for the F ex -Fl model) and (e,f) (for the 
-Pe^-Sjogren model) of Fig. [8] We can see that the be- 
haviour of the long time limit of the tracer correlator is 
in accordance with the results of Sect. IV B (see Fig. [6]). 
For the real part, the long time limit goes down to zero 
monotonously with increasing _F CX , whereas for the imag- 
inary part, the long time limit first increases from the 
zero value at F cx = and then goes down to zero, until 
the critical force value F e c x is reached (see the dashed line 
in the plots). For F ex > F£ x , the long time limit remains 
zero. 

For F cx < F c c x , the real part of the correlator is a 
monotonously decaying function of time, whereas the 
imaginary part exhibits a maximum. For F ex > i 7 ^, both 
the real and the imaginary parts of the tracer correlator 
decay faster and faster with increasing F cx and eventually 
start to oscillate. These oscillations are due to the F ex - 
dependent term w (see eqs. (115), ( 117[ ), < |120 )), which 
dominates at high forces. The oscillations thus arise for 
all bath states at larger forces (see Fig. [8j panels (a-h)). 

For the case that the bath is in the arrested state but 
the coupling between the tracer and the bath is small 
(see panels (c-d) of Fig. [8]), the long time limit of the 
tracer correlator is equal to zero for all values of F ex . 
The behaviour of 4> s (t) is similar to that for the case of 
strong probe- bath coupling for F cx > F° x . 

Finally, for the case of the liquid bath (panels (g-h) 
of Fig. [8]), we can see that besides lowering the value 
of the intermediate /3-plateau with increasing F ex (this 



effect corresponds to the effect of decreasing long time 
limit with increasing F ex for e > 0), also the time scale 
of the a-process (i. e. the final decay from the /3-plateau 
to zero) decreases. After the critical value F° x (e = 0) is 
reached, the /3-plateau becomes zero and the difference 
between the /3-process and the a-process disappears. The 
behaviour for F ex > F^ x (e — 0) is similar to that of the 
case e > (as was already mentioned in this section): 
the overall time scale of the decay decreases and for large 
enough F cx , the tracer correlator starts to oscillate. 

In the next sections, we consider the /3-relaxation and 
the a-relaxation regions of the i^^-Sjogren model in more 
detail and make some quantitative predictions. 



E. The /3-correlators 

In this section we consider the behaviour of the correla- 
tors around the /3-relaxation plateau and perform a (non- 
linear) stability analysis of the arrested/localized part of 
the correlator. Here, classical MCT has provided the 
deepest insights by deriving results like the factorization 
theorem and power-law relaxation during the so-called 
/3-process. The q-dependent factorization theorem shows 
that the dynamics on all length scales follows a single, 
time-dependent function, the so-called /3-correlator G(t). 
It depends sensitively on the separation to the MCT bi- 
furcation and introduces algebraic decay into the dynam- 
ics. We perform the non-linear stability analysis in order 
to investigate the de-localization transition at finite force 
in more detail. 



1. The fi-scaling equation 
We use the ansatz 



m = f + c(t) 



(138) 
(139) 



with the assum ptions that f s , / fulfill the long-time limit 
equation ( 121 1 and the /3-correlators G s (t), G(t) are small 



\G s (t)\, \G(t)\ « 1. 



(140) 



Using the standard steps (partial integration etc.), we 
rewrite eq. (115) in the form 



d t </, s (t) = -u><t) s (t)+m(t)-^ / dt'm(t-t')<t> s (t'). 

at Jo 

(141) 



We insert eqs. (138), (139) into (141) and obtain 
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FIG. 8: Time-dependence of the probe-particle correlators from the schematic models with increasing F ex . The left column 
shows the real parts, the right column the imaginary parts. Panels (a-f) correspond to glass states, where the bath is nonergodic; 
panels (g,h) are for a fluid state, here F£ x = -F e c x (e = 0). Panels (a-d) are for the F e x-Fl model where the bath is completely 
arrested (</> = 1), and panels (e-h) are for the f e:E -Sjogren model. Panels (c,d) are for a weakly coupled probe (v s < which 
remains mobile in glass even at vanishing force. The parameter values are v a = 2.0 (a,b), v s = 0.8 (c, d) and e = 10 4 , v 2 = 2.0, 
v s = 4.0 (e, f), e = -10" 7 , v 2 = 2.0, v s = 4.0 (g, h). 
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O t G s (t) = - W (/ S +G S (<))+« S (f s *+G s *(t))(f+G(t))- 



dt 



dt'v s (f s *+G s *(t-t>)) (f+G(t-t>)) (f s +G s (t>)) (142) 



While eq. ( 142 ) has not been solved yet for all relevant 



cases, a number of solutions exist and provide insight into 
the tracer dynamics close to derealization. 

If we retain only the terms of the order not higher than 
linear in G s , G in eq. (142), neglect the time derivative 
and make use of eq. (121), we obtain 



w 

v. 



+ G s *(t)f(i-n 
c(t) (r ~\n 2 ) = o 



(143) 




2. Factorization theorem for fluid states in the F ex -Sjogren 

model 



FIG. 9: The critical amplitudes (eqs. ( |147[ ), ( |148[ )) of the 
-Fex-Sjogren model in the fluid state («2 = 2.0, v s — 4.0, 
F e c x (e = 0) = 0.61). 



Here we want to consider the F ear Sj6gren model in the 
liquid state (e < 0) for the case that the external force is 
smaller than its critical value for e = 0. This means that 
the /3-relaxation plateau is non-zero. 

We consider eq. ( 143 1 and choose f s to be the long 
time limit of the tracer correlator for e = 0. So, f s y^0 
) can be considered as a linear equation for 
re given G(t). Expressing u/v s in terms of 
namely /f and /| and using 
we obtain 



and eq. ( 143 
G s {t) with t 

the known functions of F ex , 
the long-time limit equation (IT 2 Th 



-c s (t) r f+G s (ty r f (i-n+G(t) \r\ 2 (i-n = o 

(144) 

In terms of Gf (t), G 2 (t), i. e. the real and imaginary parts 
of G s (t), eq. (144| can be easily solved with the result 



Gf(t) 

cm 



= hiG{t), 
= h 2 G(t), 



(145) 
(146) 



of the force values F ex < F£ x . At F ex = 0, the func- 
tion h 2 (F ex ) starts linearly from zero, whereas hi(F ex ) 
starts quadratically at a non-zero value, as required by 
symmetry. 

To check our results numerically, we plot the (3- 
correlators (defined in Eqs 
and then scaled according to t 
in Fig 



138 



we 

139])), hrst 



unsealed 



10 



le expressions ( 145 1, ( 146 ) 



We see, that Gf(t), G 2 {t) indeed collapse on 
the master curve given by G(t), if one is not too far away 
from the plateau. This holds for more than ten decades 
in time in Fig. |10| 



The critical correlators 



In this section we want to consider the /3-correlators at 
the critical force value F ex — F£ x . So, we set f s = in 



where 



hi = 



eq. ( 143 ) and obtain 

G s {t)uj + G s *(t) vj = 0. 



/(4/j 



(150) 

1/1 (2/T/l +(2/r + ,/'| - f! ) (1 - A s )) ( U7)As can be readily seen, in terms of Gf (t), G%{t) the above 

:2 f? (i - 7f) - (2/f + if - ff) UP - 7nf ■ 



equation can be rewritten as 



hi f (/| 2 



/f') + l/ s | 2 (l-/f) 



2//r/i 



(148) 



J 



cm 
cm 



0. 



(151) 



This result generalizes the factorization theorem of MCT 
to forced probes. As a check, we set F ex = and obtain 
the well-known result of the Sjogren model [TBI |3"5] : 



where the matrix J is identical with the one given by 
eqs. (133), (134). Thus, except from the trivial solu- 
tion Gjtt) = G s 2 (t) = 0, the solution of eq. (fl~50| ex- 



i 



(149) 



ists only when the determinant of the matrix J vanishes, 
which is exactly the bifurcation condition, as discussed 
in Sect. IV C If it is fulfilled, the solution of eq. (150) is 



The critical amplitudes hi , hi are plotted in Fig. [9] as 
functions of F ex for fixed values of v% and v s . Both func- 
tions increase monotonically in the (meaningful) region 



not unique and represents a relationship between Gf(t) 
and G s 2 (t). 

In our case the codimension of the bifurcation, i. e. the 
dimension of the critical space is one, as was shown in 
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FIG. 10: Beta correlators, unsealed (upper panel) and scaled 
(lower panel) of the -F e:E -Sj6gren model in the fluid state for 
F cx < F c c x (e = 0) (v 2 = 2.0, v B = 4.0, e = -10" 11 , F c c x (e = 
0) = 0.61). 



Sec. IV C Thus at the bifurcation, the correlators G\(t), 
G|(i) are proportional to the critical eigenvector of the 
stability matrix and thus to each other. 

Considering also the next-to-leading (quadratic) terms 
in eq. (142) and still assuming f s = 0, we obtain 



v. 



d 

dt 



dt' G s {t')G s *{t-t'). 

(152) 

For the F ex -Fl model, we have to set / = 1 and G = 
in eq. (152 1. The left-hand side of eq.(152) vanishes 



due to the condition ( 150 ) and we obtain the following 



equation for the critical /3-correlator: 
d rt 



= 



dt 



dt' G s (t')G s *{t-t'). 



(153) 



This equation can be solved by means of the power-law 
ansatz 



G s (t) = t x + it y . 



(154) 




le+10 



FIG. 11: The critical correlators of the F ex -Fl model for dif- 
ferent values of v a . The dashed curve represents the asymp- 
totic power law (|157l|. 



We get under the integral in eq. ( 153 1 the expression 



t' x (t - t') x + t ,y (t - t') y + i [t ,y (t - t') x - t ,x (t - t'f] . 

(155) 

Using the identity 



d 
dt 



dt' (t-t') x t' y = t x+y 



T{x + l)T{y + l) 



T{x- 



1) 



(156) 



where T(x) is the gamma function, we see that the imag- 
inary part of the right-hand side in eq. ( 153 1 vanishes. 



The choice x — y = — 1/2 lets also the real part of the 
right-hand side in eq. ( 153 ) vanish, since then the denom- 



inator in eq. ( 156 ) diverges 



We thus found the power law solution 
G s (t) = t- 1 ' 2 + it- 1 ' 2 



(157) 



of the equation ( 153 ), which gives the critical /3-correlator 
of the i^x-Fl model. Fig. 11 shows the critical correlators 
for different values of the parameter v s . We see that 



the power law (157) for the real and the imaginary part 
indeed holds asymptotically for large times. The solution 
( |157 1 can still be multiplied by an arbitrary prefactor. 
This expresses the scale invariance of the eq. (1531. The 



correct prefactor can be found by matching to the initial 
decay. 

If we consider now the Fex-Sjogren model, we get from 
eqs. (p2l, ([1501 



G s *(t) G(t) 



d 
dt 



dt' G s {t')G s *{t-t'). 



(158) 



This equation has no simple power-law solution, since 
with the power-law ansatz, the imaginary part of its 
right-hand side would vanish, whereas the left-hand side 
would still have a non- vanishing imaginary part. 
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F. The a-relaxation 

We want to consider now the a-decay region of the F ex - 
Sjogren model for e < (see pa nels ( g-h) of Fig. [8] and 
the discussion at the end of Sec. IIVDI) in more detail. 



The a-relaxation behaviour of the schematic model 
without the external force is a well known example of 
the second scaling-region of MCT, describing the final 
decay of the correlator on time scale r, the so-called fi- 
nal, or a-relaxation time. For the present discussion we 
recall, that the second relaxation step of the correlators 
asymptotically (for e — !• 0~) follows a scaling- function 



<p(t) ~ j> 



t 



(159) 



As we saw in Sec. IV D the presence of the external 
force influences both the time scale of the a-decay and 
the height of the /3-plateau. So, a simple scaling law like 
(159 1 cannot work any more. However, if we rescale the 

by the factor ^wg^p 



2 decay from the same 
we see that 
and 



12) 



amplitude of the correlator tfif 
(i = 1, 2), so that both <p\ and 
(F ex — real part-) plateau (see Fig 
with increasing _F ox , the shape of both am iy 2 
varies slightly. There is also some difference between the 
shape of 4>\ and 4>2 & t the same value of F ex , which de- 
creases with increasing F ex so that for F ex = 1.4 (slightly 
below the critical force) 0f and 0| almost match. 

This observation justifies us to propose the (approxi- 
mate) generalized ansatz 



-curves 



s (*)~rcF ex )^ 



t 



r s {F ex ,e) 



(160) 



suggesting that there is still a universal decay function 
but accounting for the change in the plateau value. <p s (t) 
is suggested to be real, which means that both the real 
and the imaginary parts of <p s have the same shape. The 
precision of (1601 can be considered as acceptable if one 



realizes that the change of the decay time scale with in- 
creasing F ex by several orders of magnitude has a much 
stronger effect than the minor change in the shape of the 
curves. 

To determine the tracer a-time scale r s , we match 



(160) to the /3-decay law (see Sec. IVE2) 



4>m = f!+h i G(t) = f? 




(161) 



where the asymptotic form of the bath beta-correlator 
G(t) = — (i/r(e)) b was used [15], to obtain 



F ex )=r(e) 



1/6 



(162) 
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FIG. 12: a-decay of the correlators, rescaled in the amplitude 
according to their /3-plateau values. The parameters are v s = 
6.0 , v 2 = 2.0, e = -10 -11 (this corresponds to F£ x = 1.445). 



where i = 1, 2 corresponds to the real and imaginary 
part, respectively. This result will be used in the next 
section to analyse the low-force behavior of the tracer 
friction coefficient in a fluid host. 



G. The tracer friction coefficient 

Within the framework of the schematic models, where 
no wave vector dependence of the correlators is present, 
we define the friction coefficient increment following our 



considerations in Sec. II C as the time integral over the 
product of the real part of the tracer correlator and the 
bath correlator: 



dt<t>\{t)<l>{t) 



(163) 



1. F ex -Flmodel 

We start by considering first the F ex -F\ model, since 
exact analytical results are available here. Eq. ( 120 ) 
reads in the Laplace space: 

-i(z^(z)+l) = -(l-iF ex )^(z)-v s (z^{z)+l)^*(-z) 

(164) 

We use the following definition of the Laplace trans- 
form 

f{z)=LT[f(t)](z) = i ^ dte izt /(*), (165) 
Jo 

with the properties 

LT[d t f(t)}(z) = -i( z f(z) + f(0)), (166) 



LT 



dsf(s)g{t - s) 



LT[f*(t)} (z) 



-if(z)g(z), 
-(/)*(-*)■ 



(167) 
(168) 
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For the calculation of the friction coefficient, only the 
imaginary part of (j) s (z — 0) is of interest, since 4> s (z — 

0) = - ST + * Jo " (*)■ So > we have A C = 

= 0)}. We set z = in cq. (1641, then the 
product z <j) s vanishes and for 



K{^(z = 0)}, 
^ 2 = 3{^(z = 0)} 

we obtain the following system of equations 



which yields 



1 + W S 



> S 2 = 



-V a S l 



l + Vs 

FL + 1 - vl 



ii 



(169) 
(170) 

(171) 
(172) 



(173) 



So, we have an exact analytical result |39j and see that 
the friction coefficient exhibits thinning behaviour with 
increasing F ox . 

As expected from symmetry, AC starts out quadrati- 
cally for small external forces for v s < 1, i. e. for the case 
of low probe-bath coupling. For v s > 1, i. e. for strong 
probe-bath coupling, expr. (1731 can be rewritten as 



1 



F 2 - F< 



2 : 



with 



(174) 



(175) 



(according to eq. (131) with / = 1). Note that eq 
applies only if F ex > F£ x , otherwise the tracer is lo< 
and A( = oo holds. Due to the identity F 2 X - F£ x 2 = 
(F ex - F^ x )(F ex + F£ x ), A( diverges at F£ x according to 
the asymptotic power law 

A C ~- 



pc 



(176) 



In Fig. [13] we plot the numerical and the analytical val- 
ues of AC for different v s and observe quite a reasonable 
agreement. The deviations increase with v s (following 
the general trend in the numerics to become unstable at 
higher values of the probe-bath coupling strength) and 
can be considered as a quality measure of the numerical 
procedures used. 



2. F ex -Sj6gren model 

For the F ea -Sjo gren m odel, the results of the numerical 
integration of eq. ( 163 ) are shown in Fig. 14 (continuous 
lines) as function of F cx for different values of e. We ob- 
serve thinning with increasing force and see that for large 
forces, all the curves collapse on the same limiting curve. 




FIG. 13: Friction coefficient increment from the F ex -Fl 
model. The circles show the numerical values, and the 
continuous lines show the analytical values calculated from 
eq. fl73l). 



This limiting curve corresponds to the large F ex limit of 
the schematic model, where the w-term, which contains 
F ex dominates and the memory term can be neglected so 
that one obtains a 1/F 2 X decay law for A£. This behav- 
ior also agrees with the low-density approximation in the 
microscopic MCT equations (see eqs. (101 102)). 

For e > 0, the curves diverge at the critical value 
of force, which increases with e. On the fluid side 
for e < 0, two different decay regimes can be distin- 
guished: the strong decay from the initial (linear re- 
sponse) plateau for F ex < F^ x (e = 0) and the further 
decay for F ex > F£ x (e = 0), which approaches the e = 
limiting curve with increasing force. The initial decay for 
e — > 0~ can be analyzed analytically using the a-decay 
law derived in Sec. |IVF[ since the time integral over 
the correlators is dominated by the a-decay region for 
F ex < F£ x (s = 0). Relations ( |l63| , Jl59| , ( p0| lead to 

AC * f?(F cx ) J Q °° dtp {Tfib^)^ (4y) • From FigS 
(g) we see that (j) s {t) decays much faster than (j)(t) and 



thus a further approximation is justified, where <J)(t) is 
considered as constant under the integral. This gives the 
scaling 



A,(--.F„) x rL-)./7(F, s )( ^(F„) 



i/b 



(177) 



where the use of eq. ( | 162 ) was made. The inset in Fig. 14 
demonstrates the validity of the factorization of the e- 
and the F ex -dependence in A((e,F cx ): rescaling of the 
amplitude of the AC vs. F ox -curves for different e leads 
to their coincidence. To check the F cx -scaling, we plot 

( r \ 1/b 

the function ff (F ox ) ( (F cx ) J (circles on Fig. 
value of the bath beta-scaling exponent b = 



14 



the 



o3 was 

used) and observe a good agreement wi th th e results of 
the direct numerical integration of eq. (163). For small 
F cx , the expansion C(e, F cx ) = Co(e)(l ~zpt + 0{F^)) 
is valid with c ; 



8.35 (see inset in Fig. 14) 
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FIG. 14: Main panel: friction coefficient increment from 
the Feaj-Sjogren model (v2 = 2.0, v s — 4.0). Results from 

The bold black line 
-10" 3 , -10~ 4 , 
3 in ascending 



the numerical integration of eq. ( 163 1 
corresponds to e 



0. The values of e are 



, -10- 
Inset: 



b , -lo-', -io- s , 

data for e = — 10~ 



0. 



10"' 
-10" 



10" 



-10"' (for F ex < 



-1GT 
order 

F£ x ) scaled by r(e) to collapse onto the mastercurve eq. ( 1771 
(continuous line). The dashed line corresponds to the small 
F ex expansion (o(e = — 10~ 7 )(1 — c.F e 2 x ) (with c = 8.35). 



V. SUMMARY AND CONCLUSIONS 

The main objective of this work was to extend the stan- 
dard mode-coupling theory for the motion of a tracer par- 
ticle in a dense colloidal suspension near the glass transi- 
tion to the case, where the tracer experiences an external 
force F£ x , which cannot be assumed to be small compared 
to the internal units of the system. This means, an at- 
tempt is made to go beyond the linear response regime. 
We use formally exact generalized Green-Kubo relations 
and follow the ideas of the integration through transients 
approach |17j . recently developed for sheared systems. 

The presence of the external force leads to a dras- 
tic difference compared to the linear response case: the 
tracer density correlator becomes complex. This is the 
consequence of the fact that instead of the unperturbed 
Smoluchowski operator Qq, the full operator contain- 
ing the external force enters 0*. fit turns out to be non- 
hermitian with respect to the equilibrium average, as the 
consequence of the fact that we consider an open system. 
Interestingly, in the mode-coupling theory under shear 
the wavevector dependent density correlators can be de- 
fined so as to remain real, as the affine drift motion of 
the particles can be taken into account rigorously. In the 
present case of force driven microrheology the drift mo- 
tion results from the particle interactions and manifests 
itself in a phase factor which needs to be calculated, and 
which turns the correlator complex. 

Despite this qualitative difference in the transient 
structural relaxation, in fluid states an external force 
gives thinning behavior of the friction coefficient akin 



to shear-thinning in flow. The difference to flow-driven 
macrorheology becomes evident in the existence of a crit- 
ical force in glass states. At F c c x a continuous bifurcation 
transition of the long-time limit of the tracer density 
correlator occurs. For F ex > F£ x , the long-time limit 
becomes zero and thus the cage surrounding the tracer 
breaks. The probe-particle becomes delocalized and can 
be pulled through the suspension. The critical force in- 
creases with the bath packing fraction (Fig. [5]). Below the 
critical force, the distortion of the cage around the (still) 
localized probe can be visualized from its formfactor in 
real space (Fig. Q. 

In Sec. |IV] we constructed the arguably most sim- 
ple schematic models by considering only two wave vec- 
tors parallel to the external force. Two different models 
are considered: the "-F ea: -Sj6gren model" , extending the 
Sjogren model [35] for the tracer coupled to a bath and 
the u F ex -Fl model", extending the Fl model of stan- 
dard MCT, which was used to describe the tracer in a 
matrix of immobile particles (the Lorentz model). The 
long-time limits and the phase diagrams could be calcu- 
lated analytically. The bifurcation at F£ x was shown to 
have codimension one. 

Numerical and asymptotic solutions of the time- 
dependent equations of motion of the schematic models 
enabled predictions of the force dependence of the tracer 
friction increment A£. Generally, the thinning behaviour 
is observed similar to the shear thinning in macrorheol- 
ogy. For the F ex -Fl model, an exact analytic expression 
for A£ could be derived, showing a power law divergence 
with the exponent —1 at the critical force. For the F ex - 
Sjogren model, scaling laws could be given for small and 
large external forces. 
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Appendix A: The adjoint of the full Smoluchowski 
operator 



The backward Smoluchowski operator (see eq. ( |15[ )) of 
the closed system is known to be self-adjoint with respect 
to the equilibrium-weighted scalar product (denned in 
eq. Q): 



(A\nlB) = (nlA\B) 



(Al) 



To calculate the adjoint of the full operator fit from 
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eq. ( 14 1 we note that 



and with eqs. (A3 1, (A4) 



(A | 0*5) = (A | SV B) + (A | Afttfl) = (ft^ | + | Afit.B) 

(A2) 



holds. The term containing Afit in (A2| reads 
(A | AfitB) =D s JdT V eq A* (F ex ■ d s )B. 



(A | AtfB) = -D s J dT V eq B (F ex • (F s + fl a ) A* 



(A3) 



(A6) 



Eqs. ( |A6|[A2|) lead us to the result (A \ QtB) = ((fij - 
Using partial integration with respect to r s and the fact Ds ' + I B > and thus to the identification 

that F ea; is space-independent, we obtain 

y dr* e9 A* (F ex -d s )B = -J drB(F ei -a s )(v*) 

(A4) 

With the relation d s ^ eq = F s vJ/ e(;j we have 
(F ex ■ 8 S ) (V eq A*) = ^ eq {F EX ■ a.) A* + ^ eq {F ex ■ F s ) A* 



n°* = nl-D s F ex • (p, + ds) 



(A7) 



for the adjoint of the full backward Smoluchowski opera- 
tor with respect to the equilibrium- weighted scalar prod- 



(A5) uct. 
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